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The chirally unbroken and the superconducting 2SC and CFL phases are investigated in the chi- 
ral limit within a Dyson-Schwinger approach for the quark propagator in QCD. The hierarchy of 
Green's functions is truncated such that at vanishing density known results for the vacuum and 
at asymptotically high densities the corresponding weak-coupling expressions are recovered. The 
anomalous dimensions of the gap functions are analytically calculated. Based on the quark propaga- 
tor the phase structure is studied, and results for the gap functions, occupation numbers, coherence 
lengths and pressure differences are given and compared with the corresponding expressions in the 
weak-coupling regime. At moderate chemical potentials the quasiparticle pairing gaps are several 
times larger than the extrapolated weak-coupling results. 
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I. INTRODUCTION 



Strongly interacting matter as it existed in the early 
universe, resides in the interior of compact stellar ob- 
jects or is produced in heavy-ion collisions is subject to 
extreme conditions. An understanding of the different 
phases of such matter in terms of the fundamental de- 
grees of freedom of QCD in the realm of low tempera- 
tures and high quark densities is still missing although 
such knowledge would be of fundamental interest. At 
sufficiently high density and low temperatures strongly 
interacting matter is expected to be a color supercon- 
ductor, which has re-attracted a lot of interest in re- 
cent years (for corresponding reviews see 0, Due 
to asymptotic freedom, it can be systematically studied 



at asymptotically large densities in a weak coupling ex- 
pansion 13 . At densities that are relevant for the interior 
of neutron stars one is in the strongly coupled regime, 
and such an expansion is lacking. Investigations 0| are 
usually done within Nambu-Jona-Lasino-type models in 
mean-field approximation. The main objective of these 
studies is the identification of possible scenarios for the 
color-superconducting state. Since the results are rather 
model dependent on a quantitative level an approach, 
which is directly based on the QCD degrees of freedom 
is highly desirable. 

Furthermore, although it has been argued on general 
grounds that quarks become deconfined at high densi- 
ties, a corresponding detailed picture is missing. This re- 
lates to the fact that already at vanishing temperatures 
and densities many aspects of the confinement mecha- 
nisms remain elusive, in spite of recent progress (see e.g. 
ref. j^]). Relations between the confining field configura- 
tions and the Gribov horizon have been established 
Those are in turn related to the infrared behavior of 
QCD Green's functions, either in Coulomb gauge 7J or 
in Landau gauge [1] QCD. In the Coulomb gauge quark 
confinement results already from an effective one-gluon 



exchange picture Il0| and confinement of colored com- 
posite states such as 'diquarks' can be studied |T3|- A 
detailed understanding of the infrared part of the Yang- 
Mills sector is still lacking, however. On the other hand, 
in the Landau gauge the gluon propagator has been de- 
termined quantitatively by different methods with con- 
sistent results, see e.g. 113 [l3> llS llal and references 
therein. The infrared behavior of the vacuum gluon and 
ghost propagators, established by these and related stud- 
ies, have provided the basis for investigations of the quark 
propagator as well as extensions to finite temperature for 
the pure Yang-Mills sector ,13- The results underline 
non-trivial strong coupling of ghosts and gluons even at 
arbitrary high temperatures. 

It seems natural to suitably extend these methods 
to finite densities especially for the strongly coupled 
regime, and compare to analytical results obtained in 
the weak coupling limit. In the present work we pro- 
ceed in this direction by extending the truncation scheme 
of the pertinent Dyson-Schwinger equations (DSE) pro- 
posed in ref. ■ We implement momentum- and energy- 
dependent self energies, which encode information about 
the non-Fermi liquid behavior of the chirally unbroken 
phase at large quark densities. This is important for the 
analysis of superconducting gap functions already in the 
weakly coupled regime We consider medium effects 
on the quark-quark interaction in the random-phase ap- 
proximation similar to the Hard-Dense-Loop (HDL) ap- 
proximation and thereby recover the corresponding an- 
alytical results in the weakly coupled regime. The pre- 
sented results allow for a quantitative estimate for the 
validity of the weak coupling results at a given chemi- 
cal potential. This is the main objective of the present 
investigation. 

This paper is organized as follows: In sect. II the trun- 
cated quark Dyson-Schwinger equation at non-vanishing 
chemical potential is derived. In sect. Ill the ultraviolet 
behavior of the gap functions is extracted analytically. 
In sect. IV we provide expressions for occupation num- 
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bers, diquark coherence lengths and the effective action. 
In sect. V we discuss the relevance and application of 
Luttinger's theorem for relativistic matter. In sect. VI 
we give numerical results for the unbroken phase. In 
sect. VII results for two superconducting phases (2SC 
and CFL) are presented, including an estimate of the 
pressure from the Cornwall- Tomboulis-Jackiw (CJT) ac- 
tion. Sect. VIII closes with some concluding remarks and 
an outlook. 



where Dc^ = — ^ff^J corresponds to the charge- 
conjugate of the covariant derivative = 9^ + igA^. 
Using the charge-conjugation matrix of Dirac spinors, 
C = 7274, the 8A^c^/-dimensional bispinors ^'(a:;) are 
defined as 



(2) 



II. QCD GREEN'S FUNCTIONS 



A. The quark propagator at non-vanishing 
chemical potential 

For the description of superconducting systems, it is 
advanta geo us to work within the Nambu-Gor'kov for- 
malism |2flj. In the following we will only consider the 
chiral limit and thus treat all flavors equally. Extend- 
ing the conventions and notations of refs. |23|, the 
unrenormalized QCD Lagrangian in Euclidean space at 
finite quark chemical potential /i and vanishing quark 
masses reads 



Note that the components of these bispinor field opera- 
tors are not independent of each other |2^. In the fol- 
lowing Nc = 3 and Nf — 3 will be used. 

In a covariant gauge, the renormalized quark DSE with 
appropriate quark- wave-function and quark-gluon- vertex 
renormalization constants, Z2 and Zip, respectively, is 
then given by 



S'^^^ ^, Al] = S [V-, V^, Al] + ^74mV 
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5-1 (p) = Z2So'{p) + Zifj:{p), 



where 



(3) 



-ip- 7 — i{p4 + ifi + gA4)j4 






-ip- 7 - i{p4 -ifj.- gAl)j4 



(4) 



r 



is the inverse bare quark propagator in the presence of a 
static, isotropic and homogeneous time component A4 of 
the gluon field. The quark self energy S(p) is given by 

= - /(0ir(^LV(g)r^G,.('?>p)^;i'(9-p). (5) 

Here D'^^{k) is the gluon propagator, and the bare and 
full quark-gluon vertex in the Nambu-Gor'kov basis is 
defined as 



.(0)a _i9 ( l^X 







NG^l. 



2 V -7mA"' 



"pa 



r-(9,p) A-^((7,p) 



(6) 



with A-^,(g,p) - -CA^(-p,-g)^C and r-^,(q,p) = 
— Cr^(— p, —q)^C. Flavor indices have been suppressed, 
they will be discussed in detail below. 



In a fixed gauge, the expectation value of A4 is deter- 
mined by its DSE, i.e. by its equation of motion. In a 
covariant gauge this corresponds to a vanishing expecta- 
tion value of 

SS'^''^\ih,ib,Af,] , ^ , , , 

' = -{dJ^' + 9r''Al{x))Ft^{x) 



igr'"'d^Cb{x)ccix), 



(7) 



where c{x) and c{x) is the ghost and antighost field, 
respectively. For a static, isotropic, homogenous, even- 
parity and T-symmetric ground state \ in a fixed co- 
variant gauge one has 



{n\DfF^^ixm) 

gr'^{n\A:{x)F;^(x)\n) 



(8) 
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and, since the ghost propagator is real and symmetric, 

(d^Cb{x)cc{x)) = . (9) 
From this one eventually obtains 

p'^ix) cx Itvd,cJ,ng {S{x)r''^^^,) = . (10) 

Therefore the static gluon fields ensure color neutrality 
(p°(x) = 0) and can be implicitly determined by this 
condition [24. 125l l26l|. For the symmetries stated above, 
this is equivalent to the exact solution of the equations 
of motion. For the quantities discussed in the following 
the neutrality condition only slightly modifies the results, 
and thus it will be neglected for simplicity. 

B. Gluon propagator and quark-gluon vertex 

The quark propagator can be determined from eq. (jS)) 
as a functional of the gluon propagator and the quark- 
gluon vertex. In the Landau gauge, which will be cm- 
ployed in the following, these Green's functions have been 
investigated for the chirally broken phase at zero chemi- 
cal potential by DSE studies and by lattice calculations, 
see e.g. HI El El El HI HI IM 113 • In the Landau 
gauge, the gluon propagator is parameterized by a single 
dressing function Z{k'^), 

DfAk') - ^a,{5,.-^-^)^. (11) 

The quark-gluon vertex possesses in general twelve lin- 
early independent Dirac tensor structures (see e.g. ref. 
|23|) with the vector- type coupling (x 7^ being the only 
one present at tree level. It has been shown that some 
of the other tensor structures a re q ualitatively important 
when studying the quark DSE jHESl- For simplicity, in 
this exploratory study the quark-gluon vertex is assumed 
to be of the form 

Tl{p,q) = ^gT{{p-qf)^,^. (12) 

The unknown function r(fc^) has been chosen such that 
the quark propagator is multiplicatively renormalizable 
an d ag rees with perturbation theory in the ultravio- 
let Eg- It can also be determined from quenched lattice 
results of the quark and gluon propagator js^l or finite- 
volume results for the corresponding DSE's in combina- 
tion with lattice data for the quark propagator (33i |. Be- 
yond the vector-type coupling contained in H12|l at least 



the term prop ortional to + is of qualitative signifi- 
cance [TfilllSj. In the chiral limit such a scalar-type term 
is, however, only non-vanishing in the chirally broken 
phase because such a contribution violates chiral symme- 
try. As we are interested here in superconducting phases 
this scalar contribution will not be present, at least, in 
the chiral limit. Therefore, the use of the approximation 
(|12|l is certainly sufficient for an exploratory study. 
In the following only the product 

«.(fc') = ^'^Z{e)T{k% (13) 

will enter the quark propagator DSE. Here Zip is the 
quark-gluon coupling and Z2 the quark-wave-function 
renormalization constant, respectively. We will refer to 
as(fc^) as the effective strong running coupling because, 
especially in the framework of DSEs, it is a possible non- 
perturbative extension of the coupling into the infrared. 

It is, of course, expected that the gluon propagator 
and the quark-gluon vertex undergo changes when a non- 
vanishing chemical potential is introduced. Note that in 
superconducting phases anomalous vertices, the AJ^ in 
eq. ©, are induced. These are neglected in the study 
reported in this paper. Thus, taking into account the 
matrix structure implied by the Nambu-Gor'kov formal- 
ism and the approximation H12|) we will use 

naMD - r((p-9)')ri^L;. (14) 

with r(fc^) taken from quenched vacuum studies. 

Since we are primarily interested in chirally unbro- 
ken phases at non-vanishing chemical potential, it is im- 
portant to incorporate medium effects like damping and 
screening by particle- hole excitations. In this paper we 
add the in-medium polarization tensor with 'bare' quark 
propagators to the inverse gluon propagator. This is not 
self-consistent and needs to be investigated in further 
studies. Nevertheless, the resulting quark DSE turns out 
to be a generalization of the HDL approximation, with 
the important difference that the infrared behavior of 
the gluon propagator and quark-gluon vertex is nontriv- 
ial. Furthermore, even in the superconducting phases, 
the assumption of bare quark propagators will a poste- 
riori turn out to be much better suited than employing 
the vacuum quark propagator, i.e. the one reflecting dy- 
namical chiral symmetry breaking. 



The renormalized medium polarization tensor is generically given by 



-Q-mcd ab 



Up) 



\ZlF 



(27r) 



rTr 



D,c,f,NG 



q)S{p-q) 



(15) 
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For the first part we employ the approximation of vanishing quark self-energies and the vertex H14() . Then this 
expression can be straightforwardly reduced to 



_^., 2.^,a.(p^) |_0_Tr,,([7,5o(<z)>^o(p-9)]-[...],.o)- (16) 



With help of projectors transverse and longitudinal to 
the medium 



-I /I. 



0, 



PiPj 



and 



respectively, this polarization tensor can be expressed by 
two functions, 0(1^1,^4) and F{\j)\,p4), 



ah TjT 



cab T-)L 



(17) 



The evaluation of this functions is well known pertur- 
batively, and in the present case only the coupling is 
replaced by the running coupling. For small external 
momenta, the result is 



Fi\p\,Pi) = 



\P\ 



1 



«P4 
\P\ 



f ipA 


iPi 


\\P\J 


\P\. 



-o 



2m\p'Y-^ 

_ iPi^ ( iP4 
\P\ ^ \ \P\ 



with Q{x) = ^ In and 

m2(p2) = Nfa,ip^)tJ.yTr. 



(18) 



(19) 



(20) 



Adding the medium polarization to the inverse gluon 
propagator leads to 



/ T~\vac 



tab 



p 



ip') 

2 



mcd ab 



Hp))'' 



,2 



r + Gip) 

2 



P, 



\ Z{p' 

p^ + F{p)^ f"") p2 



(21) 



Taking into account this form of the in-medium polariza- 
tion, Debye screening and Landau damping are included, 
similar as in the HDL approximation. As it is phrased 
here, it becomes evident that these phenomena have a 
non-perturbative origin and require the knowledge of the 
infrared behavior of Schwinger functions, in particular of 



C. The truncated quark DSE at non-vanishing 
chemical potential 

Summarizing the above considerations we arrive at a 
truncated, self-consistent DSE for the quark propagator. 
Treating the components in the Nambu-Gor'kov basis 
separately we define 2] 



'^"^^^ " I S^ip) 



Sip) = 
E(p) = 



S+{p) T-{p) 

T+ip) S-ip) 

S+(p) $-(p) 

cf>+(p) E-(rt 



with 



A-{p) = -C{A+{-p)fC 
for e {S^,S^,T^}, 



(22) 



(23) 



S+(p) + I]-(p) = 74 (S+(p) + S-(J3))S4 and 

<t>'{p)=j4{^+{p))h4 (24) 

for a real action. One then obtains [s^ 

-1 



1 ^ ^±-1 

-1 



— Z2SQ + ZipT,^ 



± 



(25) 



combined with 

E+ip) = 



^+ip) 



ZiF^ J (2^)4 

^ a.(P)Pj, «.(fc2)P,i \ 
k'^ + G{k) B + F{k) j ' 

ZiF J (27r)* 

F + G{k) fc2 + P(fc) ] ■ 



(26) 
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Figure 1: Two versions of the strong running coupling in the 
numerical solution of the DSE (I26II . The full line represents 
the coupling ai(k^) extracted from the coupled gluon, ghost 
and quark DSE's [l^ and the dashed line aii{k^) extracted 
from quenched lattice data together with an abelian approx- 
imation for the quark-gluon vertex js^ . 



Note that the running coupling as(fc^) also enters via the 
functions G(|fc|,fc4) and F{\k\,k4). On the other hand, 
this running couphng is the only input into eqs. (|26|) . 
In the following we will present results for two different 
couplings. Both couplings are derived in a MOM reg- 
ularization scheme. One will be the running coupling 
determined from the coupled gluon, ghost and quark 
DSEs jl8j, denoted by ai{k'^) in the following, where 



a/(fc2) = 



ln(e + ai {k^/Aff' +bi {k^Aj) 



(27) 



with a/(0) = 2.972, oi = 1.106, 03 = 2.324, 61 = 0.004, 
&2 = 3.169 and A/ = 0.714GeV. As the quark-gluon ver- 
tex used in this study contains a sizeable "scalar" con- 
tribution (which is due to dynamical chiral symmetry 
breaking) the corresponding running coupling together 
with the approximation (|12|l underestimates the chiral 
condensate significantly but should be, as mentioned 
above, the appropriate choice in chirally restored phases. 
We consider this coupling as a lower bound. On the other 
hand, the running coupling extracted from quenched lat- 



tice data together with an abelian approximation for the 
quark-gluon vertex serves as an upper bound and 
will be denoted by a//(fc^) in the following. As the cor- 
responding parameterization is quite lengthy, we refer to 
eqs. (8)-(16) of ref. 32J for its definition. Both couplings 
are shown in Fig. ^ 

In order to obtain a self-consistent solution of the sys- 
tem of equations (|26|l with H19|) . H2UI) and the running 
coupling as input, we will first consider the color and fla- 
vor structure of the quark self energy. We will restrict 
ourselves to scalar pairing and treat the 2-flavor color- 
superconducting (2SC) and color-flavor locked (CFL) 
phases separately. In the chiral limit, the pairing struc- 
ture can be described by a symmetric matrix M such 
that the projectors on its eigenspaces Pi together with 
the matrices Mi :— MPi form a closed basis under the 
transformations Pi XaPiKi and Mi X^MiXa, re- 
spectively. Therefore, we can parameterize the r.h.s. 
of eqs. H26|l by the renormalization point independent, 
Dirac-algebra valued self energies Y.f (p) and (fi^ (p) : 



% 



(28) 
(29) 



For an even-parity and T-symmetric phase those are 
given by [s^ 



S+(p) = -iplY.\ip)-i^^Y.lip) 



(30) 



= 74 E K^^pW?^ 

e=±l 

^tkv) = (74^01,(P)+0J,,(P))75 

= 75 E (31) 



e=±l 



where we made use of the positive and negative energy 
projectors Ai = \ (\±i-i4j), ■§ ^ p ■ p = 
tpp = bjp^A and LUp — ipi + ^. 

Introducing positive parameters 5i > via the de- 
composition M^M = J2i ^iPi the quark particle-particle 
propagator can be written as 



Z2S+ip) 



(-ZP4 - m) (1 + S+ ,(p)*) + e\p\ (1 + SX .(P)*) 



\{ipA - m)(i 



KM) 



e\p\{l 



K,{pW+^^\<l^eAp)\' 



(32) 



where we have made use of the relations TiF,i{p) — Y,F,i{~p) ■ The zero of the numerator at p4 = defines 
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the quasiparticle Fermi momenta pp , and the zero in the 
denominator provides the corresponding dispersion rela- 
tion. To first approximation, the energy gap in the 
excitation spectrum is therefore given by 



A? 



ip) 



(33) 



\p\=Pp ,P4=Q 



III. 



ULTRAVIOLET FINITENESS OF THE GAP 
FUNCTIONS 



Similar to the ultraviolet analysis of the quark mass 
function in the chirally broken phase [s^j we determine 
here the ultraviolet behavior of the gap functions. For 
large external momenta p, such that p^, /^^ ^ 

p^ and m? <gC p^, the breaking of Lorentz covariance is 
negligibly small, and thus the self-energies ^"^{p) and 
are to a very good approximation functions of the 
four-momentum squared, p"^, only. Furthermore, in the 
denominators of the integral kernels the self-energies can 
be safely neglected. Also the medium modifications of 
the gluon propagator will then be insignificant, and the 
gap equations reduce to 



(34) 



^ ^ kl + {q-kf\\ 



X p-g 1-2 



fc2 



(35) 



As the running coupling is a slowly varying function for 
large momenta it is safe to apply the angular approxi- 
mation Q;s(fc2) « as{p^ 0{p'^ - q^) + e{q^ - p^)). The 
remaining angular integrations can then be done analyt- 
ically. The gap function ■ {p) decreases for large p^ at 
least like times logarithmic corrections, and thus to 
order one has 4>\i[p) ~ 0. The equation for the 

gap function 0^ ^ [p) then reads 



16n 



dq 



^ as{q^)<P+,{q) 



(36) 



The color-antitriplet channel is attractive and for this 
channel we get X^MiXa = — f-^i as compared to AaA^ — 
^ = ACp- The color sextet channel is repulsive. For this 
one we have A^M^ Ao = The anomalous dimensions 

of the gap functions 7^^^, to one- loop order can then be 
read off from the coefficient in eq. (|36|l . Comparing to the 
corresponding anomalous dimension of the mass function 
7m = 12/(33 — 27V/), they are given by 



and 



7^,3= 7m/2= 6/(33 -2A^/) 



70,6 = -7m/4 = -3/(33 -2A^/) 



(37) 



(38) 



in attractive and repulsive channels, respectively. Similar 
to the chiral quark condensate, see e.g. ref. [38j, one is 
now able to define a renormalization-group independent 
diquark condensate. The asymptotic behavior of (f)^ ^{p) 
is given by the so-called 'regular form' 



1 



In 



A2 



(39) 



IV. RELEVANT QUANTITIES 

A. Occupation numbers and diquark correlations 

Once the quark propagator is known, one can extract 
number densities, occupation numbers and the diquark 
coherence lengths. Within the Euclidean formalism, the 
number density p is calculated as the derivative of the 
generating functional of the connected Green's functions 
with respect to the chemical potential /i. For the homo- 
geneous phases, considered here it is given by 

p - Z2{{i^Hx)i,{x))) 

= lim -Tr^i.cj^ATG (^2 74 lArG'5(x4,0)) 

d^'p f dpi ^^^^^ ^ ^ ^^^^ 



E 



(27r)3 
9i 



(27r)3 



27r 

d^P ni{p), 



(40) 



where gi — 2rank(Pi) is a degeneracy factor, and the 
occupation numbers ni{p) read 



nt{p) 



Z2 
An 



dpi TiD (74S'/(p4,p)) 



(41) 



The p4-integration has to performed first to make this 
expression well-defined. The relation between the density 
and the Fermi momentum will be further discussed in 
sect. El 



7 



The quark-quark correlation lengths provides a mea- 
sure of the size of the paired diquarks. They can be 
determined from the anomalous propagator 

T+ix-y) = mxfCM,i;{y))) 



(42) 



For a given pairing pattern, M,;, the coherence length 
is defined as 



Jd^p\T+{pW ■ 



(43) 



B. The effective action 

Studying different phases, the question naturally arises 
which phase is energetically preferred. To study this, we 
estimate the corresponding pressure difference by em- 
ploying the Cornwall- Jackiw-Tomboulis (CJT) formal- 
ism ^39j] . which provides the effective action F as a func- 
tional of the expectation values of fields and propagators 
in presence of local and bilocal source terms. In par- 
ticular, for QCD in the Nambu-Gor'kov formalism the 
functional dependence on the quark propagator is given 



T[S] 



Tr„.D_c,/,iVGLn5 ^ 



1 

■^'TT:p,D,c,f,NG (l - Z2Sg^S) + T2[S] 



(44) 



this theorem is based on the fact that the functional F2 [S] 
is invariant under shifts in the momentum, i.e., 



ST2[S] — —-ZipTrpD^cj.NG 



d ' 



= 0. (46) 



Using Z274 «) Ing = igfx 



(5^1 - ^ifS), eq. ^ and 
the spectral representation of the propagator one can 
show that 



1 



d^p 
(2^ 



P4=o- 



27r 



-Tr 



D,cJ,NG 



Log(S-'] 



(47) 



P4=0+ , 



where the cut in the complex logarithm has been, as 
usual, put onto the negative real half-axis. This, for a 
single fermion species, amounts to 



— TrzjLog (5-1) 



P4=0- 



Tr^Log(5-i) 



Pi=0+ 



^ In (iLOpC + \/IVTs2) 



P4=Q- 

P4=0 + 



2 if D{p,p4, = 0) < 
if D{p,p4 = 0) > ' 



(48) 



where D = + w^C^ -I- B^. Since the Fermi surface 

in a Fermi liquid is defined by D{p = pp,p4 — 0) = 0, 
it is natural to extend this definition to a sign change in 
D{p,p4 — 0). For a gapped mode this corresponds to a 
singularity. Finally, we note that 



Tr^GLog(5-i) = Log {S+ 



(49) 



Here r2[iS] is the 2-particle irreducible part of this ef- 
fective action. Note that the self energy can be ex- 
pressed as the functional derivative with respect to r2 [S] , 
ZipJ^[S] = —26r2[S]/SS, and thus possesses a corre- 
sponding functional dependence on the full quark prop- 
agator. One can "integrate" the DSE and obtains an 
approximate, yet thermodynamically consistent, effective 
action: 



L2[S] ~ "^Trp,i5,c,/,AfG(l 

+ const. . 



Z2Sa'S) 



(45) 



as weh as D{p,p4 = 0) of S+ and {Z2S^ + Zip^-) 
change sign at the same momenta, since S*^ is only in- 
vertible iff Z2{Sq)~^ + ZipT,^ is. In an isotropic phase 
one therefore concludes that 



(50) 



Especially, the Fermi momenta p^.i are calculated by the 
zeros and poles in deto^cj (5"'") |p4=o- We can therefore 
easily determine the density as long as the symmetry 
^fAp) = ^F,ii—p)* in our numerical solutions is not 
spontaneously broken. 



V. ON LUTTINGER'S THEOREM 

Luttinger's theorem can be summarized as follows: 
Provided the fermion propagator is positive at the Fermi 
energy, p^ = 0, the volume of the Fermi surface at fixed 
density is independent of the interaction. The proof of 



VI. RESULTS FOR THE UNBROKEN PHASE 

Before turning to the superconducting phases and in 
order to provide a case for comparison, the unbroken 
phase in the chiral limit is studied first. In eqs. H2t)|) 
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the functions and accordingly T^, are set to zero. 
This yields a self-consistent integral equation for the self 
energy S"*", which we solve numerically. 

To understand the behavior of this numerical solution, 
we derive an approximate form for its imaginary part. 
First, the inverse propagator is decomposed into positive 
and negative energy parts 

S+{p)-' - 74At5+(p)-i+74AZ5^(p)-i . (51) 

The quasiparticle propagator 



Z2' -Wi 1 - 



ImI]I|l(p) 
Pi 



-Ai-|p1+ReS+(p) 



(52) 



near the Fermi surface has been subject of several weak 
coupling analyses [111 . Especially the imaginary part 
of the self energy, encoding the wave-function renormal- 
ization on the Fermi surface has been studied in detail. 
Within the present truncation one is able to general- 
ize the investigation described in ref. such that the 
momentum dependence of the wave-function renormal- 
ization is taken into account. For bare quarks, corre- 
sponding to a 1-loop approximation, the contribution of 
the transversal gluons to the self energy is given with 
I = \p \ ~PF by 



Kip) 



4i 
'3^ 



dkA I dk arctan 



P4 + ^4 + il 



(53) 



The integrand is discontinuous for |?| < A:, and the leading 
non-analytic contribution can be extracted from 



.ImS^(p) 



P4 



^KiP) 



d 
dpA 
Ai 



(54) 



dk- 



as{k'^)k 



k>\l\ P -I- fTO2(p)lM- 



The main contribution to the integral comes from scales 

1/3 

of the order k ~ (tt m^(fc^)|p4|) . Note that fc ^ for 
P4 — > 0. Employing this implicitly given scale one can 
approximate 



ImS+(p) 4a,(fc2) 



Pi 



9 



In 



PF 



^m^e)\p,\ 



A3 



(55) 



where Ajjv is an ultraviolet cutoff. This demonstrates 
the well-known fact that the long-range (static) gluon 



interaction renders quark matter in the unbroken phase 
into a non-Fermi liquid. In the present context it is obvi- 
ous that this non-trivial feature depends on the infrared 
behavior of the product of gluon propagator and quark- 
gluon vertex. 

Here also a qualitative difference between the couplings 
displayed in Fig. ^ becomes important. Since the cou- 
pling aiiq^) becomes almost constant for q < 0.5GeV, we 
can estimate that k becomes negligibly small, i.e. fc ~ 0, 
if 



\Pi\ < 



(0.5GeV)3 lOMcV 



For the coupling a//(q^) this is not possible because the 
coupling vanishes in the infrared. One would need \p4\ <C 
(10MeV)/(^2[Qev2]) to study this effect, which is far 
beyond the scope of the present work. 

In Fig. 121 we compare the non-analytic p4 dependence 
at the Fermi surface to the corresponding numerical re- 
sult, employing the coupling ai{k'^). One clearly sees, 
that if k and Auv are chosen accordingly, the approx- 
imation H55|) works well. In addition, we display the 
dependence on \p\ (which is neglected in weak coupling 
analyses) demonstrating that the singularity only occurs 
at IpI = pp. 

For completeness, we present numerical results for the 
self-energy functions and in the unbroken phase, 
see Fig. 121 Displayed are their values at the Fermi surface 
as a function of the chemical potential fj, and the momen- 
tum dependence for p4 = at = 0.5GeV and lOGeV. 
One sees that, in the unbroken phase, the approximation 
using bare quarks in the medium polarization is justified 
a posteriori since, as expected, EJ^, <^ 1. 



VII. RESULTS FOR 2SC AND CFL PHASE IN 
THE CHIRAL LIMIT 

In this section we combine the presentation of the nu- 
merical results for the scalar 2SC and CFL phases for 
three massless flavors. 

The color-flavor structure of the 2SC phase is given by 
M — T2 which leads to 



Z2 

z 



IF 



■ (l>2SC ^2®T2. 



(56) 



Here (r, g, b) and (w, d, s) label color and flavor. The ma- 
trices Ai and Ti are the Gell-Mann matrices of color and 

flavor space, respectively, completed by Aq = tq = -y/fl- 

As described above, Pi, denote projectors on the coupled 
color-flavor space. 

The CFL phase is given by the ansatz M — X2® T2 + 
A5 (81 T5 -I- A7 (8) T7. Note that M is antisymmetric in color 
and flavor, respectively. The corresponding condensates 
can be chosen to be invariant under a3(8)3 = 1®8 
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Figure 2: The logarithmic singularity on the Fermi surface 
as a function of p4, compared to the approximation given in 
eq. 1)55^ at scales k = and k — fi, (upper panel) and its mo- 
mentum dependence for fixed P4 ^ (lower panel) . Here the 
coupling ai{k^) is employed. Auv ~ l.SGeV is chosen such 
that the analytical approximation fits the numerical result. 



transformation with generators Tq — j4^]. Then the 
matrices {Pi} turn out to project onto the irreducible 
representations. We therefore obtain 




(57) 
(58) 



S={0,1:3,4,6,7,8 



'TS 



where we have also introduced the commonly used 'an- 
titriplet' and 'sextet' pairing function 0^" = -I- |(/)g 

and (j)^ ^ + i^s- 

Both phases have been studied in a weak coupUng anal- 
yses in the HDL approximation. Including in this approx- 
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Figure 3: The self energies T,\ and in the unbroken phase 
at the Fermi surface as a function of the chemical potential 
fi (upper panel) and as a function of the momentum |p| for 
jj, = 0.5GeV and lOGeV for pi — (lower panel). Note that 
for p4 = all dressing functions are real. 



imation the "normal" quark self energies, the quasiparti- 
cle gap at the Fermi surface is given by [l9| 



weak A 



= 512 tt" 



1 i = 2SC 

2-1/3 i^i 



(59) 



The momentum dependence takes a more complicated 
form but for g'^ In (^jJ- / (t>^g^ak i) O'^^ '^^'^ neglect the 
quark self energies and obtains 



(IpI) 



(60) 



cos 



■In 



weak.i 



3V2Tr \\p-f^\ + 4iP) 



For the weak-coupling approximations of the gap func- 
tions displayed in the following figures we use the iden- 
tical running coupling as the one employed to obtain 
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Figure 4: The quasiparticle gaps 4>'^ and anti-quasiparticle 
gaps at the Fermi surface for the 2SC (upper panel) and 
CFL (lower panel) phases for the coupling ai{k^). These are 
compared to the extrapolated weak coupling result 4't)eak- 
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Figure 5: The quasiparticle gaps 4>'^ and anti-quasiparticle 
gaps 0" at the Fermi surface for the 2SC (upper panel) and 
CFL (lower panel) phases for coupling aii{k^). These are 
compared to the extrapolated weak coupling result (jiweak- 



the displayed numerical solution, respectively. Note, 
however, that in the strong coupling regime the weak- 
coupling expressions ceases to be valid. 

The results for the gap functions at the Fermi surface 
the for 2SC and CFL phases are shown in Fig. 0] for the 
coupling a/(fc^). In Fig. [S] the corresponding results are 
shown for the coupling aii(k^). In both figures it is evi- 
dent that the extrapolated weak-coupling results for the 
value of the quasiparticle gap differs by a factor of more 
than five for fj, « SOOMeV, resulting in quasiparticle- 
pairing gaps larger than 50MeV. It is worth noting that 
the results are less sensitive on the coupling than the 
mass function in the chirally broken phase. This can be 
attributed to the fact that a stronger coupling also re- 
sults in stronger screening and damping. In addition, 
also the anti-quasiparticle gap function are determined 
self-consistently. Nevertheless, it is remarkable that the 
ratio of the 2SC and anti-triplet CFL gaps is almost equal 
to 2^^^ as in the weak coupling analysis. The reason for 
this behavior is that the normal self energies are only 
weakly modified compared to the self energies in the un- 



broken phase. This amounts to an effective "decoupling" 
of normal self energies and gap functions, and thus the 
factor 2^/"^ is directly inherited from the pairing pattern 
in color- flavor space. 

As can be inferred from Fig. |31 and the estimate given 
in eq. the energy gap in the excitation spectrum 

is at most 15% smaller than the gap function at the 
corresponding kinematical point. (Note that the corre- 
sponding coefficients 6i in the decomposition of M^M 
are given by: \/S2SC = — 1 and \/Ss = 2.) We also 
emphasize, that 4i^{p) — > 4>^ {p) for ^ 0, which can 
be non-zero, indicating a Bose-Einstein condensation of 
diquarks. This is found for the coupling a//(fc^). How- 
ever, this is not expected to be the energetically favored 
state in the vacuum, since spontaneous chiral symmetry 
breaking has not been taken into account here. 

In Fig. El we present the momentum dependence of the 
quasiparticle gaps (/)"*" and anti-quasiparticle gaps 4>~ for 
P4 = for several chemical potentials, obtained with the 
coupling Q;/(fc^). For large values of /i the gap function is 
concentrated around the Fermi surface and the quasipar- 
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Figure 6: The momentum dependence of the quasiparticle 
gaps (j)^ (upper panel) and anti-quasiparticle gaps <j)^ (lower 
panel) for pi = Q normalized to their value at the Fermi sur- 
face. 



tide gaps show a cosine-like behavior near the maximum. 
However, at chemical potentials of the order Aqcd, the 
Fermi surface is no longer the main contributing region 
to the gap integral, and even the maximum of the gap 
function is no longer on or close to the Fermi surface. 

As described in sect. IIV Al the occupation numbers of 
the quasiparticles are calculated from the normal propa- 
gators. These are displayed in the upper panel of Fig. [7| 
for the 2SC phase, employing the coupling Q!/(fc^). The 
occupation number of the gapped red and green u- and d- 
quarks is, as expected and due to the pairing, a smooth 
function around the Fermi momentum. On the other 
hand, the occupation number of the blue u- and d-quarks 
changes rapidly. Within the present approximation, i.e. 
neglecting quark self energies in the medium polarization, 
these quarks are in a state at the borderline between a 
Fermi and a non-Fermi liquid. 

The decoupled s-quarks are in the unbroken phase and 
one clearly sees, especially when comparing to the blue 
u- and d-quarks, that the occupation number changes 
smoothly. As already mentioned, this is another feature 
of non-Fermi liquids which, strictly speaking, indicates 
the breakdown of the quasiparticle picture. We also find 
a significant depletion at |p| = due to the interaction. 
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Figure 7: The occupation numbers (upper panel) and the 
function defined in eq. l(iH . pL,i{v), (lower panel) for the 
different quark channels in the 2SC phase at jj, = IGeV with 
the coupling aj{k^). As can be seen in the panels all Fermi 
momenta are (almost) equal and given by pF,i ~ .937GeV, 
c.f. the dashed-dotted line in the lower panel. The lines of 
the different channels are nearly indistinguishable. 



The density can either be obtained by integration of 
the occupation numbers or via the Luttinger theorem 
(see sect. IV|l . The conditions for the latter are fulfilled 
since all self-energies are real for p4 = 0. To analyze this 
remarkable feature we define the function 



dqq^ni{q) 



(61) 



which has to obey the limiting behavior 
lim^^oo PL.i(f^) — PF,i- From Fig. [7| one sees that 
this limit is assumed for v > ppi, and that the expecta- 
tion derived from Luttinger 's theorem is nicely fulfilled 
within the numerical accuracy Note also that all 

Fermi momenta are very close to each other, a fact which 
justifies a posteriori the approximation of neglecting the 
neutrality conditions. 

It is instructive to compare the coherence length of the 
diquarks, ^-)-, to the mean-free path as determined from 
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Figure 8: The ratio between the coherence length and mean- 
free path in 2SC and CFL phase for the coupling ai(fc^) (up- 
per panel) and aii{k^) (lower panel). 



the density. The results are shown in Fig. |S1 as a func- 
tion of the chemical potential. Although the two different 
couplings lead to a distinctive pattern for these ratios it 
is safe to conclude that the size of a Cooper pair at mod- 
erate chemical potentials, /i « 500MeV, is only several 
times the mean-free path, the precise value depending on 
the diquark channel and the employed coupling (similar 
results were already presented in ref. (44] for the 2SC 
phase.). Although there is an analogy to the crossover 
between a BCS-type superconductor in weak coupling to 
the strongly coupled regime (Bose-Einstein condensate) 
the size of the diquarks suggest that the 2SC and CFL 
phases resemble a strongly coupled BCS system. How- 
ever the result indicates, that mean-field approximations 
are very questionable. In the limit /x — s- 0, we again find a 
different behavior for the two employed couplings indicat- 
ing Bose-Einstein condensation of diquarks for a//(fc^) 
for very small chemical potentials. 

Finally, the pressure difference between the supercon- 
ducting phases (2SC and CFL) and the unbroken phase 
are displayed in Fig. |51 This quantity is given by the 
negative effective action difference, and is here calculated 



> 
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Figure 9; The pressure difference between the superconduct- 
ing phases (2SC and CFL) and the unbroken phase employing 

within the approximate C JT formula given in eq. H44|) . In 
the chiral limit, the CFL phase is the most favored phase, 
as expected. This is in good qualitative agreement with 
the parameterized form 

Ap = J2ranm)^^0-. (62) 

i 

which is (slightly) a generalized formula as compared to 
the one obtained in the weak coupling regime [40|. 

VIII. CONCLUSIONS AND OUTLOOK 

We have presented the solution of a truncated Dyson- 
Schwinger equation for the quark propagator at non- 
vanishing chemical potential in the chiral limit. The 
truncation includes the in-medium modification of the 
gluon propagator and has been chosen such that, at 
vanishing density, known results for the vacuum and at 
asymptotically high densities the corresponding weak- 
coupling expressions are recovered. We have determined 
analytically the ultraviolet behavior, i.e. the anomalous 
dimensions, of the gap functions. 

The most remarkable feature of these self-consistent 
solutions are the values of the gap functions: The 
quasiparticle-pairing gaps are several times larger than 
the extrapolated weak-coupling results even for sizeable 
chemical potentials /i ~ IGeV. This is true irrespective 
of the considered phase (2SC or CFL) and the employed 
running coupling. 

The investigation, presented here, provides the start- 
ing point for more realistic studies involving finite current 
masses. Using different quark masses the question of im- 
plementing neutrality conditions (which could be ignored 
here) becomes important. 

Due to solutions of the gluon- and ghost Dyson- 
Schwinger equations at non-vanishing temperatures [T^ 
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an investigation of different regions of the phase diagram 
opens up. However, to solve the complete coupled system 
of the Dyson-Schwinger equations for the QCD propa- 
gators remains a challenge. A study of the question of 
deconfinement due to a non-vanishing quark density re- 
quires certainly also a better understanding of the in- 
frared behavior of these Green's functions and thus an 
infrared analysis of the coupled Dyson-Schwinger system 
at non- vanishing temperatures and densities. 
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